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Abstract 


This report summarizes the advances made to date in near-wall Reynolds-stress closures. 
All closures examined are based on some form of high-Reynolds-number models. Researchers 
argue that because turbulent diffusion is small near the wall, its modeling is relatively insignificant 
compared to viscous dissipation and velocity-pressure-gradient correlation in the near-wall region. 
Consequently, most near-wall closures proposed to date attempt to modify the high-Reynolds- 
number models for the dissipation function and the pressure redistribution term so that the resultant 
models are applicable all the way to the wall. Furthermore, the near-wall closures examined solve 
a turbulent kinetic energy dissipation-rate equation to complete the closure. Some solve the 
equation that governs the transport of the dissipation rate, while others solve an equation that 
governs the transport of the actual dissipation rate minus its wall value. Again, the high-Reynolds- 
number form of the equation is modified to predict the near-wall flow. The near-wall closures are 
examined for their asymptotic behavior so that they can be compared with the proper near-wall 
behavior of the exact Reynolds-stress equations. A comparison of the closures' performance in the 
calculation of a low-Reynolds-number plane channel flow is carried out. In addition, the closures 
are evaluated for their ability to predict the turbulence statistics and the limiting behavior of the 
structure parameters compared to direct simulation data. It is found that three near- wall Reynolds- 
stress closures give the best correlations with simulated turbulence statistics; however, their 
predictions of the near-wall Reynolds-stress budgets are incorrect. A proposed modification to the 
dissipation-rate equation remedies part of those predictions. Therefore, further improvements are 
required if a complete replication of all the turbulent properties and Reynolds-stress budgets by a 
statistical model of turbulence is desirable. 
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1 . Introduction 


Statistical models of turbulence have made significant advances since the mixing-length 
model of Prandtl 1 . This early model assumes gradient transport and prescribes, on the basis of 
empirical data, an algebraic behavior for the mixing length. The corresponding velocity scale is 
given by the product of the mixing length and the mean shear rate. Subsequent models improve on 
the empirical input and propose to solve either a turbulent length scale or a turbulent velocity scale 
equation instead of relying on a prescribed mixing length. Further improvements can be achieved 
by solving the length and velocity scale equations simultaneously with the mean flow equations. 
Consequently, the length or the velocity scale does not have to be tied to the mean shear rate. In 
other words, small scale turbulence is not required to have a direct link to the large scale motion. 
The length and velocity scale equations can be derived rigorously from the one-point statistical 
equations. However, in order to close these equations, the flow Reynolds number is assumed 
large. This assumption renders the models invalid for near-wall flows and, therefore, limits the 
range of applicability of the models. As a result, these so called two-equation closures are not very 
suitable for turbulent flows with such complexities as adverse pressure gradient, streamline 
curvature, fluid rotation, recirculating and reattachment behavior, etc 2 . 

Efforts to improve these closures have been carried out in two directions. One is to extend 
the two-equation closures to low-Reynolds-number flows 3 . Another is to relax the gradient 
transport assumption and to solve the one-point Reynolds-stress equations directly 4 . This latter 
approach gives rise to the commonly known Reynolds-stress or second-order closures. Again, the 
Reynolds-stress equations are closed by invoking the high-Reynolds-number assumption. As 
such, they are not valid for near-wall flows. These two improvements give mixed results 2 . 
However, they still fail to capture the important characteristics of near-wall turbulence 4 and the 
closures are performing poorly for such complex flows as turbulent curved-pipe flows-’A In 
order to further improve the performance of second-order closures, a number of researchers have 


proposed to extend the second-order closures to near-wall flows 7 ' 15 . In general, the extension is 
based on well established high-Reynolds-number models for the various terms in the Reynolds- 
stress equations. A common approach is to modify the pressure redistribution model in order to 
account for the "echo" effect in the proximity of a wall 7 - 9 * 10 - 12 . Some of these second-order near- 
wall closures have been tested for a variety of simple and complex turbulent flows, and they show 
improvements for the flow cases that have been calculated 4 * 7 ' 23 . However, the asymptotic near- 
wall behavior and budgets of the modelled Reynolds-stress equations have not been properly 
analyzed. As a result, the internal consistency of the modelled equations in the near- wall region is 
in doubt and needs to be analyzed and compared with the behavior and budgets of the exact 
Reynolds-stress equations. 

With the advent of supercomputers, it is now possible to solve the Navier-Stokes equations 
directly to simulate turbulent flows 24 ' 27 . As such, the statistical approach together with its 
inherently empirical closure assumptions becomes less attractive for turbulent flow calculations. 
The direct simulation approach yields turbulence results near a wall that are in excellent agreement 
with measurements and, for the first time, offers fluctuating pressure information that is otherwise 
not available from experiments 28 - 29 . This approach provides a powerful tool for the study of 
turbulent flows. At present, it is limited by flow geometries, flow Reynolds number and the speed 
and capacity of supercomputers. Therefore, it cannot be used to tackle such practical problems as 
internal flows through turbomachines. On the other hand, the turbulence information created by 
direct simulation of simple turbulent flows would be invaluable to turbulence modellers whose 
closures could be improved by incorporating this new knowledge, in particular, the pressure 
information, into their closure assumptions. In view of the fact that direct simulation of turbulent 
flows is far from being able to compete on an equal footing with turbulence modeling in terms of 
practical applications, it would seem prudent, at least for the near term, to make use of the 
simulation data to evaluate the internal consistency of second-order near-wall closures. There are 
nine such closures 7 ' 15 at present Most of these closures have not been properly analyzed for their 


2 



near-wall asymptotic behavior and their consistency with simulation results. This suggest that, 
before further validations and improvements of these closures are carried out, a thorough review of 
their ability to predict near-wall flows is beneficial. 

The present objective is to review the second-order near-wall turbulence closures proposed 
to date. A detailed analysis of the asymptotic near-wall behavior of these closures and their 
comparison to simulated and measured behavior is carried out. Therefore, the relative strengths 
and weaknesses of the closures could be investigated and modifications proposed. It is hoped that, 
through this review, a physically and mathematically more consistent closure could be identified 
and possible improvements suggested. 
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2 . Basic High-Reynolds-Number Models 


This paper proposes to review the advances made in second-order near-wall closures. An 
examination of the near- wall closures 7 ' 15 put forward to date reveal that, in one way or another, 
they represent extension of high-Reynolds-number closures to near-wall flows. The extensions 
come in different forms. However, they all involve modifying the transport equation for the 
dissipation rate of the turbulent kinetic energy, e, and the models for viscous dissipation and 
pressure redistribution. In view of this, it would facilitate discussion of near-wall closures by first 
summarizing the high-Reynolds-number models used in these closures. The closure by Shih and 
Lumley 7 is not reviewed here, because the inner boundary conditions of their closure have to be 
applied at a streamline in the inertial sublayer and not at the wall. Therefore, it is different from 
other near-wall closures. In these other closures 8 " 15 , the integration of the governing equations 
can start right from the wall and, hence, the wall boundary conditions for all turbulence quantities 
could be satisfied exactly. The eight near-wall closures reviewed are those proposed by Hanjalic 
and Launder 8 , Prud’homme and Elghobashi 9 , Kebede, Launder and Younis 10 , So and Yoo 11 , 
Shima 12 , Launder and Tselepidakis 13 , Launder and Shima 14 , and Lai and So 15 . For ease of 
reference later, these closures are denoted by the following abbreviations; namely, HL, PE, KLY, 
SY, SH, LT, LSH, and LS, respectively. 


The governing equations for an incompressible, stationary turbulent flow can be concisely 
written in Cartesian tensor as 



= 0 , 


( 1 ) 


9U; 

u ^= 


1 3P 


p 3xj + 8x k 


9U; 

^ - UlU * 


( 2 ) 


where uj and Uj are the i th components of the fluctuating and mean velocity, P is mean static 
pressure, x k is the k th component of the Cartesian coordinates, p is fluid density and v is fluid 
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kinematic viscosity. These equations are not closed because of the presence of the Reynolds-stress 


term, -UjU k . A second-order closure of Eqs. (1) and (2) involves solving additional transport 
equations for the Reynolds stresses and e. These equations can be symbolically written as: 


Cy - DX + dJ + Pjj + cpP + Ojj - Ejj 


( 3 ) 


C e = D e +D J + P e + V-^e + ^ • 


( 4 ) 


The terms from left to right in (3) represent the convection, molecular diffusion, turbulent 
diffusion, production by mean shear, pressure diffusion, pressure redistribution, and viscous 
dissipation of UjUj . Similarly, the terms in (4) represent convection, molecular diffusion, 
turbulent diffusion, production, extra production by mean shear, and viscous dissipation of e. The 
last term £ is usually neglected while the terms Dy, Oy, ey, D £ , P e and Z) E are modelled so that 
Eqs. (1) - (4) are closed. The linkage to the wall is then provided by wall function 
approximations, such as those proposed by Hanjalic and Launder 30 and Launder et al.^ 1 

All near- wall closures put forward to date except S Y use a common high-Reynolds-number 
e-equation. The models proposed for D £ , P e and D e are given by 9[C ec j(k/e) u k Uj 3e/3xj]/9x k , 
C El (e/k)P andC e 2 £ 2 /k, respectively, where C E( j, C El and C £ 2 are model constants, P = P^/2 and k 
= UjUj /2 is the turbulent kinetic energy. On the other hand, SY solve the transport of e = e - 
2v(3k 1 ^ 2 /3x2> 2 rather than e, where X 2 is normal to the wall. Thus defined, e goes to zero at the 
wall. The models they assumed for D £ , P e and D e are 3[(v t /a e ) 3e/3x k ]/9x k , C El (e/k)P and 
C e2 e 2 /k, respectively, where v t = C^k 2 /? is the eddy viscosity and ct e and are model constants. 
In this form, the equation is identical to the high-Reynolds-number form used by Chien32. 

The model used for ey is the isotropic model of Kolmogorov 33 (hereafter denoted by KV) 
while the models adopted for Oy vary from the retum-to-isotropy model of Rotta 34 (hereafter 
denoted by RA) to the mean-strain model of Launder et al. 31 (hereafter denoted by LRR). As for 
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Dy, the models proposed range from the isotropic model of Shir 33 (hereafter denoted by SR) to 
Daly and Harlow's model 36 (hereafter denoted by DH) to Hanjatic and Launder's model 30 
(hereafter denoted by HL72). A more complicated model for Dy has also been proposed by 
Cormack et al. 37 However, a later comparison of these four Dy models by Amano and Goel 38 has 
shown that HL72 gives the most reasonable results among the four, even though the best result is 
given by solving a transport equation for the triple velocity correlation. 

The different proposals adopted for Dy, <J>y and £y are listed in Table 1. In this tabulation, 
LRR1 is used to denote the LRR model without the wall pressure "echo" term while LRR2 denotes 
the LRR model with the "echo" term included. The modified LRR model is a different version of 
LRR1 with modifications made in both the turbulent and the mean-strain part. Lumley's 
proposal 3 ^ to modify the turbulent part so that it goes to zero at the wall is adopted. This is 
accomplished by introducing a "flatness" parameter that varies from zero at a wall or an interface to 
one for isotropic turbulence. For completness sake, the various high-Reynolds-number models are 
summarized below. The turbulent diffusion models are given by: 



1(2 3 urn; 

t V 


(5a) 


DH = 


3 xl 


k d U iU j 

c < - war 


HL72 = 


3x k 


i (/ 3u;u; 3u k u; 

C s *■ I U k U£ -57-^ + UjUj -537— + UiUfi 


3xn 


3xo 


dUjU k ' 

dx ft . 


(5b) 

(5c) 


The dissipation model is adopted from Kolmogorov 33 , or 
KV= £ . 


( 6 ) 


Finally, the different models proposed for Oy are: 
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e f 2 

RA = -C^ u iUj - 36 ^ , 


LRR1 = RA + O ij>2 = RA - o^Py - Pj- Pi ^D y - 1 5jj fj- J x k Sy 

’4— , 3— 3 lf C 3 k3/2> 

LRR2 — LRR1 + ^ u k u m n k n m^ij " 2 u ^ u i n k n j " 2 U ^ U J n ^ n i J I £ 


• r 3 3 1 

+ [<Dkm,2 n k n m 8 ij • 2°ik,2 n k n j * 2°jk,2 n k n iJ 


C 3 k3/2 


CjkM 


modified LRR = C 1 (AA 2 )'/ 2 e ag+ 1.1 [a ik a kj - jSjjAj - ajje - 0.6 (Py - 58 ^ 


+ 0.6ea#W 2 ™S kt - ^ 


-0.6 [(A 2 (Pjj - Djj) + 3 a m ja n j (^mn " ®mn)l 


In these equations 


9U; DU; 

P ij = ‘L UiUk ^k + UjUk ^k 


D ij = ^W^ + 

_ m au i i 
s u -Laxj + a Xi J * 


.... m .2 5 .. 

a u"L k 3 0 ‘JJ ’ 


f 2 = 1 - exp [-C 3 Uj y/v] , 
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A 2 = a ij a ij * 

A3 = ay a^ ay , 

A = 1 - 9 (A 2 - A 3 )/8 , 

t T I 

nj = (0, 1, 0) is a unit normal vector and Cj, Cj, C 2 , C 3 , C 3 , C s , aj. Pi and are model 
constants. The three constants cti, Pj and Yi are not independent, rather they are related to a 
constant C 2 . The exact relations are given in Ref. 31. Different researchers adopt different values 
for the model constants. Therefore, the readers are advised to refer to the original papers for their 
chosen values. It should be pointed out that e should be replaced by e in all the above models 
when they are used in SY. 

It is obvious from this compilation that the high-Reynolds-number versions of these near- 

wall closures are very similar. They all used KV to model and essentially different versions of 

T 

LRR to model Oy. The difference in the modeling of Dy is slight and, as pointed out by Launder 
et al. 31 , is insignificant. As a result, all these closures can be expected to perform well far away 
from a wall for a wide variety of turbulent flows because they are basically similar to the second- 
order high-Reynolds-number closure put forward by Launder et al. 31 
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3 . Near-Wall Behavior of the Reynolds Stresses 


Sufficient data are now available to establish the following relations for the Reynolds 
stresses in the near-wall region. These data are obtained from direct simulation of a plane channel 
flow 25 and from two-dimensional channel 40 ' 45 , pipe 46 * 47 and flat plate boundary layer 
experiments 48 . If u+ ms , vj ms , w* ms , k + and -uv* are used to denote (u^/u 2 ) 172 , (^/u 2 ) 1/2 , (w 2 
yu 2 ) 1 / 2 , k/u 2 and -"uv/u 2 , then the Reynolds stresses, assuming incompressible flow, in the near- 

XT X 

wall region can be expressed in terms of y + as: 


+ 

11 „ 
rms 

= a u y + + b u y +2 + ... , 

(8a) 

+ 

V rms 

= a v y +2 + b v y +3 + ... , 

(8b) 

W rms 

= a w y + + b w y > 

(8c) 

k + = 

a k y +2 + b k y +3 + ... , 

(8d) 

-uv + 

= a uvy +3 **■ b uv y ••• » 

(8c) 


where u^ is the wall friction velocity, u, v and w are the fluctuating velocity components along the 

stream (x), the normal (y) and the transverse (z) direction, the a's and b's are constants for fully- 

developed turbulent flows, and y + = yu^/v. The leading coefficients of (8) can be determined from 

existing data 25,40-45 by examining the behavior of the Reynolds stresses in the region, 0 < y + < 

10. Whenever possible, mean lines through the data are used to approximate the behavior near the 

wall. The results thus obtained are tabulated in Table 2 together with the maximum values of u rms , 

v + w + k + and -uv* (denoted by a subscript m) and their respective locations (denoted by 
rms’ rms’ 

( y +) m ). The Reynolds number Re is based on centerline velocity and channel width (or pipe 
diameter). 
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From this tabulation, it can be seen that there is a lot of scatter between the various 
measurements. For example, a u varies from a low of 0.24 to a high of 0.45 with an average value 
of ~ 0.35. It is now known that the measurements of Eckelmann 4 ^ and Kreplin and Eckelmann 42 
are not accurate in the near- wall region because they have not corrected their hot-wire outputs for 
wall proximity effects. The measurements of Kasagi et al. 49 are similar to those of Ref. 29 
because their hot-film outputs are also not corrected for wall influences. Furthermore, the 
measurements of Sirkar and Hanratty 50 give an a k = 0.05, while the data collected by Derksen and 
Azad 51 give a range of 0.025 to 0.05 for a k . These values are substantially lowered than those 
reported by Kim et al. 25 and Nishino and Kasagi 45 . Some of the discrepancies are due to 
difference in Re. However, according to the analysis of Nishino and Kasagi 45 , it seems that most 
hot-wire or hot-film measurements are affected by various error sources in the near-wall region. 
Some errors, like heat conduction and natural convection effects, could be corrected but others are 
unaccountable. As a result, the hot-wire or hot-film data are not as accurate as the particle tracking 
measurements of Ref. 45. Since the data of Ref. 45 are in close agreement with the simulation 
result 25 , it would seem that these two sets of data are most suitable for validating second-order 
near-wall closures. 

The dissipation rate of k is defined as 


e ~ v ll)ii|] ■ (?) 

If the expansions 

u = aj y + a 2 y 2 + ... , (10a) 

v = b! y + b 2 y 2 + ... , (10b) 

w = ci y + c 2 y 2 + ... , (10c) 
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where aj(x,z,t), bj(x, z, t) and Cj(x, z, t) are random functions, are assumed in the near-wall 
region, then continuity requires bj = 0. A simple manipulation of (10) then gives a k = 
v 2 ^ + cf]/2i4 and b k = v 3 ^a 1 a 2 + qc 2 where the overbar is used to denote time average. 
Substituting (10) into (9) gives the near-wall behavior of e and the result is 

e + = 2a k + 4b k y + + ... , 0 *) 

where e + = ev/u^. Therefore, the limiting wall value of e is e + = 2a k and its slope at the wall is 
given by 4b k . In view of this, it is important to determine the coefficients a k and b k correctly. 
Nishino and Kasagi 45 give a b k = - 8.4 x 10' 3 . They have also determined the b coefficients for 
the other Reynolds stresses. For completness sake, these values are quoted here as b u = 
- 0.21, b v = -5.03 x 10 -4 , b w = -0.010 and b uv = - 4.9 x 10' 5 . The negative b k value indicates that 
-E decreases away from the wall. In other words, -£ reaches a maximum at the wall. This trend is 
consistent with direct simulation results 2 **. In view of the small value of bj^, it is rather inaccurate 
to try to determine from the simulation data^5. Nevertheless, an attempt has been made and the 
value so determined is b k = -6.6 x 10‘ 3 , which is in fair agreement with that reported in Ref. 45. 
Therefore, the slope of e + at the wall is 4b k = -0.0264 and is not so small as to be negligible. 
Consequently, the assumption that 8£ + /5y + = 0 at the wall is incorrect and cannot be used as a 
boundary condition for £ as suggested in Ref. 3. The correct boundary condition should be £+ = 
2a k . In general, a k is not known a prior. However, since dVk + /dy + = Va k , this suggests that an 
alternative boundary condition could be E + = 2(5Vk + /9y + )2, a condition first derived by Jones and 
Launder 52 . Similarly, k + /y +2 = a k , therefore, £ + = 2k + /y +2 could also be used as a boundary 
condition. This means that k + /y +2 e + is exactly 0.5 and could be used to assess the closures’ ability 
to replicate near- wall flows. The b k value thus determined plus the values for the coefficient "a" 
listed in Table 2 and the k + /y +2 £ + value are used to assess the internal consistency of the eight 
second-order near-wall closures. 
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Finally, in the near-wall region, the dissipation function Ey and the Reynolds stresses UjUj 
have to asymptote correctly to the following ratios 53 ; namely. 


e ll 

_ e 33 _ 

e i3 e 

“ k ’ 

(12a) 

u 2 

w 2 

uw 


£21 

e 12 

<N 1^ 

ii 

(12b) 

vw 

uv 



e 22 

II 


(12c) 


Furthermore, Ey has to contract to 2e, otherwise (9) will not be satisfied. Here, 1 denotes the 
stream or x-direction, 2 and 3, the normal and transverse directions, respectively. Together, (8), 
(1 1) and (12) define the correct asymptotic behavior for u*uj , k, e and Ey in the near- wall region. 

The validity of any second-order near-wall closure should be judged by its ability to mimic these 
characteristics, especially when its basic high-Reynolds-number model has been extensively 
verified as in the eight cases listed in Table 1. 
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4. Near-Wall Reynolds-Stress Closures 


With this understanding of the near-wall behavior of the Reynolds stresses, the 

improvements made by various researchers to model near-wall turbulence can now be discussed. 

v 

It should be pointed out that since the terms and do not need modeling they are included in 

all near-wall modelled equations. In other words, (3) and (4) are solved without neglecting and 

D V in these equations. With the exception of LSH and LS, all near- wall closures examined neglect 
E t 

<!>”, or, it can be said that pressure diffusion effects are being modelled together with Dy by such 

models as SR, DH and HL72. Recognizing the importance of pressure diffusion in the near-wall 

region 28 * 29 , LSH and LS attempt to incorporate this effect into the closure by considering the 

modeling of the velocity-pressure-gradient correlation rather than pressure redistribution alone in 

the near-wall region. Consequently, the models that need modifications are those associated with 

D ti’ ’ °'j * and e ij‘ 

HL Closure? 

HL argue that, at least to the lowest order, v does not appear in Oy. Therefore, the model 
for Oy need not be modified for near-wall viscous effects. They also argue that, for y + < 15, the 
effects of diffusive transport on the stress budgets are insignificant. In other words, HL72 can be 
used without modification. This only leaves £y for modification. They argue that as the Reynolds 
number approaches zero, the energy-containing and dissipation range of motions overlap and £y 

could be approximated by Ey = ujuj (e/k)^ 4 . Consequently, they propose to modify £y by 


£ij - 3 £ 


U 1 U L 


[(1 - f S )8jj + ( 2k y 3 ) f s j 


(13) 


where f s = (1 + Rj/lO)' 1 and R T = k 2 /ve. This £y model contracts to 2e everywhere in the flow. 
They further modify their £-equation by defining the term Z) £ as Q 2 f e £e/k, where f £ = 1 - (2/9) 
exp[-(R T /6) 2 ]. By introducing £ to replace one of the £ in £) e , the term approaches zero at the wall 
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by virtue of ( 8 ) and (11). In addition, they neglect the effect of £ and model the extra production 
term in the e-equation to account for mean-strain generation of e, or 


\jr = C e 3 v-UjU k 



(14) 


PE Closure? 

T 

PE propose to modify the high-Reynolds-number model for by multiplying it by f^ = 
exp [-3.4 (1 + Rt/ 50)' 2 ]. They adopt a form for e,j very similar to (13), but with minor changes. 
The rationale is to provide better representation of the anisotropy of the normal stresses. Their 
model is given by 


^ij = 3 ^ij ® k u i u j 0 ^ij) f. 


(15) 


where the summation convention does not apply. This modification does not allow to contract 
to 2e near a wall. To further account for normal stress anisotropy near a wall, they propose to add 
a term, 

^ij.w = c 3 ( p ij * D ij) exp {-C 4 k 1 / 2 y/v], (16) 

to <X>jj. Consequently, their near-wall model for <hy is given by (LRR1 + w ). The models they 

used to modify the e-equation are similar to those of HL with two exceptions. In their case, the 
extra term \y is given by 



and the high-Reynolds-number model for Z) e is modified by fy. 

KLY Closure 10 

This closure is similar to that of HL, with two exceptions. One is the high-Reynolds- 
number model adopted for Oy. They use LRR2 which include the wall pressure "echo" term. 
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Another is the model proposed for Ejj. They point out that (13) does not satisfy the kinematic 
conditions (12) and proceed to propose 

Ejj = 1 8jj e( 1 - f s ) + | f s F [uiUj + UiU k n k nj + uju k n k nj + 5y u k u ft n k n J , (18) 

to replace (13). This modification satisfies (12), but it does not contract properly to 2 e if F is taken 
to be unity. However, its proper contraction is assured if F = (1 + 2.5 v 2 /k) _1 is assumed. Here, 
f s is defined differently and is given by f s = exp [-Rj/40]. The E-equation used is identical to that 

of HL. 

SY Closure 11 

SY essentially follow the arguements of HL and suggest the modification of Ejj only. Their 
concern is to make the governing equations balance, at least to the lowest order of y, near a wall. 
As a result, they adopt Chien's low-Reynolds-number E-equation 32 which neglects V and proposes 

to model £ by 

£ = - exp [-C 4 yu t A'l 

r 

to account for near-wall viscous effects. Chien also defines D £ = d[(v t /o e ) (3£/9x k )]/3x k , where 
v t = C^(k 2 /£)f 2 . With this definition for v t , the model requires e to behave like y 2 near a wall, 
instead of y as suggested by (1 1). Nevertheless, this behavior renders the models for <Djj to vanish 
at the wall and allows SY to modify Ejj only. Their proposal is 


2 t ~ 2v 5j£ 5jn 

e ij = 3^ij e + y2 


One drawback of this model is that it does not contract properly to 2 e. However, the closure 
equations do have the advantage of remaining in balance as a wall is approached. 
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SH Closure 12 

SH follows the arguements of Lumley 54 to rearrange and £y. Instead of modeling 
these two terms separately, the treatment re-writes them into 


°ij = ~ 3 6 ij e + 


O 

_ 1J ’ 


2 S 

+ 


where Oy = Oy,i + <X>y ,2 and <J>y,i represents the turbulent part and d>y ,2 the mean-strain part of 
<t>y. This way, SH does not have to model £y. He proposes to model the square bracketed term in 
(21) by the RA model except that Cj is replaced by C^l- (1 - 1/Ci)f w ], where f w = exp 
[-(0.0 15k 1/2 y/v) 4 ] is specified. The term <J>y, 2 is given by (7b). In addition, SH proposes a near- 
wall correction to (21) given by 

Oij.w-f* (WPy - 1 5ijP) - O.OB^Dij - 1 5yP) + O.OSkSij] . (22) 

The £-equation is modelled by first analysing its near-wall behavior using (8d) and (1 1). This 
results in a term of 

[" 2 + 9°^] k " 2k fw ^ 23 - ) 


to effect balance of the £-equation in the near-wall region, after replacing C e 2 £ 2 /k by C e 2 eiT/k. 
Here, T is defined by SH to be T = £ - v8 2 k/3xj3xj. As for \| f, SH lumps it together with P E and 
proposes 

P e + ¥ = C el (l+f w )|P . (24) 

In view of the definition of T, the boundary condition for £ becomes £ = v8 2 k/3xj3xj at the wall. 

It should be pointed out that even though Shima 12 proposed to model d>y, w with three 
terms on the right hand side of (22), he actually neglected the second term within the square 
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bracket in his closure calculations. His justification was that it is relatively small compared to the 
first term. In view of this, the SH model used in the present calculations is the final model adopted 
by Shima; that is, (22) is implemented without the second term in the square bracket. 

LT Closure 13 

This closure adopts (18) for Ey to account for near-wall viscous effects and an identical e- 
equation to that proposed by KLY. Therefore, the only difference between LT and KLY is in the 
modeling of Oy . KLY use (7c) for Oy while LT use (7d). 

LSH Closure 14 

This closure adopts an approach similar to that of SH and proposes only to modify Oy. 
The modifications implement the suggestion of Lumley-^, so that the resultant Oy vanishes at the 
wall. Instead of just varying Cj to accomplish this objective, LSH propose to modify the 
constants Cj, C 2 , Cj and C 2 in LRR2 by 


Cj = 1 + 2.58 AAf [1-exp {-(0.0067R T ) 2 }] , 

(25a) 

C 2 = 0.75A 1/2 , 

(25b) 

c'i = -| c i + L 67 , 

(25c) 

c 2 = max [(| c 2- ^/C 2 , 0] . 

(25d) 


As for the e-equation, ^ is assumed to be zero and y is taken to be given by V - (¥1 + ^ 2 )^’ 
where 


(p \ 

V,=2.5A - - 1 , 

(26a) 

U J 


¥ 2 = 0.3(1 - 0.3A 2 ) exp[ - (0.002R t ) 2 ] . 

(26b) 
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LS Closure JS 

The rationale behind the derivation of this closure is essentially similar to that used in SY, 
namely, the first and foremost concern is to effect a balance of the modelled equations in the near- 
wall region. While SY carry this out to 0(y 2 ) for the 22 component, LS carry out the balance to 
0(y 3 ) because the direct simulation results 25 * 26 show that pressure diffusion along the 22 
component is important near a wall. In (3), the term (OP + Ojj) vanishes at the wall for all 
components while the term Ojj does not. Consequently, LS propose to model (OP + Oy) in such a 
way that the model approaches LRR1 far away from the wall. This is accomplished by adding a 
term Oy W to LRR1, where 

fw.l. (27) 


^j.w = Ci UiUj - fSijk)- | (uju k n k ni + UjU k n k n}) + a* ^Pjj - ^ Sjj Pj 


and f w | = exp [-(R T /150) 2 ]. Furthermore, Ey is modelled in a manner similar to KLY and the 
resultant model is given by 

Ey = |e (1 - f w l ) 5jj + f w l | [ujUj + ujuT n k nj + uju^ n k nj + njnj u k u 4 n k nj/ 


1 + 


3 u k u n k n ^ 
2k 


(28) 


The e-equation is modified similarly to that of SH so that de/dt is required to have the proper 
behavior at the wall. LS propose to model (y + C E j Pe/k) by C e 2 (1 + a f W 2 )eP/k, where a 
depends on Re as a result of the direct simulation study of Mansour et al. 55 The expression 
suggested by LS for a is a = [1.0 - O.bexpC-Re/lO 4 )]. Here, f w> 2 is given by f w 2 = exp[- 
(Rj/64) 2 ]. Furthermore, LS propose the following expression for namely. 


\ = fw,2 


( 7 5 C £2 -2)V r -^ 


(29) 
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where e* = e - 2vk/y 2 . This way, the wall boundary condition for e is 2v(3Vk/3y) 2 and the e- 
equation still asymptote correctly at the wall as suggested by SH. 
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5 . Near-Wall Asymptotic Behavior of the Closures 


In order to analyse the near- wall asymptotic behavior of each closure, expansions (10) are 
assumed for uj. Similar expansions can also be assumed for Uj. From the definitions of k and e, 
the use of (10) leads to 

k = |A iy 2 + B iy 3 + E iy 4 + ... , (30a) 


e = v Aj + 4vBjy + vFj y 2 + 


(30b) 


where Aj = (af + c]), Bi = (aia2 + Ci C 2 ), E } = aj + 2aja3 + bj + c^ + 2cicJ and F^ besides 
being functions of the statistics of aj and Cj, is also functions of the derivatives of a.\ and Cj. 
Furthermore, if (10) and the expansions for U, are substituted into (3), it can be easily shown that 
the lowest order term among the different components of Qj, Dy and Pjj is 0(y 3 ). Therefore, to 
0(y 2 ), they do not affect the balance of (3) in the near-wall region. The analyses of the terms D v , 

U 

(d>]j + Ojj) and Ey show that the lowest order term in (0^ + ^ 22 ) has to be 0(y 2 ) in order to 
balance the term (D^ - Ey) to this order, and the lowest order term in the other components of (Oy 


ij 


+ O.j) has to be O(y) in order to balance the same components of (Dy + Ejj) to O(y). This means 
j ij J 

that all the terms in (3) vanish at the wall. Therefore, the modelled form of (3) should also possess 
this basic property. 


Using expansions (10) and (30), the various turbulent diffusion models (5) can be analysed 
for their near-wall behavior. With the exception of SY, the asymptotic behavior of k/e is given by 
(y 2 /2v) + 0(y 3 ) and k 2 /e = Ajy 4 /4v + 0(y 5 ). Consequently, the lowest order term among the 
components in DH and HL72 is 0(y 6 ) and in SR is 0(y 5 ). As for SY, e is solved instead of £. 
Also, e is required to behave like y 2 near a wall 32 . This leads to a behavior of 0(y°) for k/e and 
0(y 2 ) for k 2 /e . Therefore, the lowest order term among the components of SR is 0(y 2 ). The 
corresponding behavior for DH and HL72 is 0(y 4 ). In other words, the asymptotic behavior of all 
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diffusion models is not correct compared to the original behavior of D-. Nevertheless, with the 
exception of SR used in S Y, the lowest order term among the different components of the various 
diffusion models is 0(y 4 ), which is smaller than 0(y 3 ) for dI. This means that in the asymptotic 
analysis of the modelled form of (3), only the term D- and the models for (O? + Ojj) and Ey need 
be considered. The result of this analysis is tabulated in Table 3 together with the behavior of the 
exact term given by (D v + OP- + Oj: - e,;). In this tabulation, the lowest order term in each 
component of the modelled and exact form of (D^ + O? + <J>y - Ey) is listed. 

From this tabulation, it can be seen that S Y, SH, LSH and LS satisfy the basic requirement 
that the modelled Reynolds-stress equations vanish at the wall. This is accomplished without 
imposing any constraints on the behavior of the random functions aj, bj and q. PE also satisfies 
the requirement for all components except 13 if Cj is taken to be unity, a condition correctly 
pointed out by Launder and Shima 14 . However, Cj is chosen to be 1.17 in PE. It should be 
pointed out that the 1 1, 33, 13 and 22 components of SH and LSH vanish at the wall as a result of 
setting Ci = 1 as suggested by Launder and Shima^ 4 . As for HL and KLY, the 11, 33, 13 and 22 
components do not vanish at the wall even when Cj = 1 is assumed. This is also true of the 13 
component of PE. The 1 1, 33, 13 and 22 components of LT are also finite at the wall irrespective 
of the value of Q. However, these expressions are too cumbersome to be listed explicitly in Table 
3. It should be noted that all the finite expressions listed in Table 3 involve the time-averaged 
correlations of the random functions aj, b, and Cj. Since the equations have to remain balance at the 
wall, certain constraints would have to be imposed by the HL, PE, KLY and LT models on the 
random functions aj, bj and cj so that the 11, 33, 13 and 22 components would go to zero at the 
wall. Therefore, this means that the near-wall turbulence statistics would be required to behave in 
a certain manner that might or might not be correct. Furthermore, according to the asymptotic 
analysis of (3), the lowest order term of (D v + 0>? + <Dy - Ey) for the 1 1, 33 and 13 components is 
0(y 3 ), while that for the 12 and 23 components is CKy 4 ) and that for the 22 component is 0(y 5 ) 
(see Table 3). In this analysis, the d/dt term in (3) is neglected because of the stationarity 
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assumption. The closures, SY, SH, LSH and LS, give an incorrect asymptotic behavior for all six 
components. However, LS tries to match the 12, 23 and 22 components to a higher order of y in 
the hope that the resultant closure could predict more correctly the anisotropic behavior of the 
normal stresses in the near-wall region. 

The failure of the closures to make the 11, 33, 13 and 22 components of (D v + <X>P. + O:; - 

ij »J J 

£y) vanish at the wall without imposing contraints on the model constants and/or the random 

functions could be traced to the incorrect modeling of <X>y or Ey. Since departures from local 

anisotropy are indistinguishable from contributions to the process of Oy or Ey, it follows that either 

or both processes could be modelled in such a way that the resultant term (D v + <J>P. + O:: - £;;) 

ij iJ J ,J 

vanishes at the wall for all i's and j's. After all, the original term does vanish at the wall for all 

components. This arguement is used in SY to accomplish the desire objective, however, at the 

expense of solving a pseudo dissipation-rate equation. In HL, Ey is not modified to compensate 

for the incorrect asymptotic behavior of LRR1. Rather, it is modified to approximate the behavior 

of £jj as Reynolds number goes to zero. Consequently, the incorrect asymptotic behavior of LRR1 

remains in the closure. The same can also be said of PE and KLY. Even though these closures 

add a Oy w term to their Oy model, the rationale is to account for the "echo" effect of pressure 

fluctuations near a wall. No consideration has been given to partially compensate for the neglected 

pressure diffusion effect and, thus, to reproduce the correct asymptotic behavior for (D v + OP. + 

ij >J 

Oy - Ey). SH and LSH recognize this shortcoming and the suggestion of Lumley 39 is applied to 
remedy the incorrect asymptotic behavior of LRR1. Even then, the behavior of the 12 and 23 
components is correct to O(y) only while the 22 component is correct to 0(y 2 ). This behavior is 
not quite consistent with the original behavior of the UjUj equations. In LS, a special attempt has 
been made to cancel out the incorrect asymptotic behavior of LRR1 near a wall. The result is a 
Oy w proposal that partially models OP near a wall, while at the same time compensates for the 
non- zero wall components of LRR1, so that the resultant modelled components of (<J>P + Oy) 
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vanish at the wall. This way, the behavior of the 12 and 23 components is improve to 0(y 2 ) and 
that of the 22 component to 0(y 3 ). 


All the dissipation-rate equations used approach their respective boundary conditions 
correctly at the wall. However, not all equations satisfy the additional wall condition found to be 
important by Shima 12 . This condition is the coincidence of 3[vd 2 k/9xj3xj]/dt and de/dt at the wall 


and it leads to 


de 

3t 


a 

at 


a 2 k 

V 

axjdxj. 


-2vFi + 24 v 2 E i 


(31) 


In the exact equations for k and e, this coincidence condition is ensured by many additional 
correlation terms. Some of these terms are neglected in modeling, therefore, it is possible that the 
resultant k and e are bound to each other in ways that are not found in the exact equations. 
Shima 12 proposes to partially compensate the effects of the neglected correlation terms by 
modeling the e-equation so that the coincidence condition at the wall is satisfied. The result is an 
additional £ term in the modelled e-equation. This idea is adopted by SH and LS. The impact of 
this on the modelled results will be assessed in the next section when the channel flow calculations 
are discussed. 
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6 . Comparisons with Plane Channel Flow Data 


The near- wall closures examined above are used to close the mean flow equations (1) and 
(2) written in Cartesian coordinates for a fully-developed channel flow at a Re = 6500. Since the 
flow is fully-developed, the mean flow equations can be reduced to a single equation involving the 
mean velocity U and the Reynolds shear stress 15 . Thus, the problem is simplified to solving one 
mean flow, four Reynolds-stress and one dissipation-rate equations. These are ordinary 
differential equations, therefore, they can be solved by a known numerical method, such as the 
Newtonian iteration technique used in Refs. 11 and 15. This means that a common, well 
established numerical technique can be used to perform all calculations and the numerical errors 
involved will be the same for each model calculation. In other words, a true comparison of the 
performance of the near- wall closures can be made for this particular problem. Any discrepancies 
in the results could be attributed to model differences and not to numerical errors. 

A fixed grid is maintained for all calculations. The grid distribution is such that there are 
five grid points in the region, 0 < y + <, 5, 15 grid points in the region, 15 < y + ^ 65, and 31 grid 
points in the region, 65 < y + < Re x , where Re x = u^H/v is the Reynolds number based on the 
friction velocity and 2H is the channel width. Since Re x is related to Re by Re x = (u x /2U 0 )Re, 
where U 0 is the centerline velocity, it serves as the only input paramenter to the problem subject to 
the following boundary conditions for U, the Reynolds stresses and e. If the x-coordinate is 
aligned with the channel mid-plane, the boundary conditions at y = H can be written as, 

U = u 2 = v 2 = w 2 = uv = 0, (32a) 

e = 2v(dVk/dy) 2 , (32b) 

while the boundary conditions at y = 0 are given by, 

de/dy = du 2 /dy = dv 2 /dy = dw 2 /dy = 0 and uv = 0. (33) 
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These boundary conditions are applicable for all closures except SY. In that case, (32a) and (33) 
apply but (32b) does not The latter condition should be replaced by e = 0, because the e-equation 
solved in the SY closure is a modified form of Chien's equation 32 , which is a transport equation 
for a pseudo dissipation rate. 

Direct numerical simulation of the Navier-Stokes equations has also been used to solve the 
channel flow problem at the same Re 25 ’ 26 . Turbulence statistics of the calculated flow are reported 
in Ref. 25 while the Reynolds-stress and dissipation-rate budgets are given in Ref. 26. As a 
result, the near-wall asymptotic behavior of the Reynolds stresses can be determined from these 
results and their limiting values are found to be consistent with another set of experimental data 45 , 
even though the two sets of data are obtained at different Re. In view of this, these two sets of 
data 25 >26,45 ^ mos t appropriate for assessing the internal consistency of the limiting behavior of 
near-wall closures. Specifically, the modelled limiting behavior of the Reynolds stresses, such as 
that given in equation (8), is compared with the chosen data. In addition, from (8) and (1 1), it can 
be deduced that k+/e+y +2 = 0.5 and (aj + a* + a* )/a k = 2.0. The ability of the near-wall closures 
to recover these limiting values is analysed. Finally, comparisons of the distributions of the 
Reynolds stresses across the channel and their near-wall budgets will also be presented. It is 
hoped that, through this rigorous comparison, the strengths and weaknesses of the various near- 
wall closures could be identified. 

The simulated, measured and modelled limiting values at the wall are shown in Table 4. 
Since Nishino and Kasagi 45 did not report on any e measurement, the limiting value k+/e+y +2 
could not be estimated. Furthermore, their value of (aj + a ^ + a^)/a k is larger than 2.0. This 
means that their measurement of a u is probably too large. Other than that, their measurements are 
very consistent with those of Ref. 25. It should be pointed out that only the modelled results of 
HL, PE, KLY, SY, SH, LSH and LS are given. In the SY calculation, HL72 and LRR1 are 
chosen as the models for the diffusion and pressure redistribution terms, respectively. The same 
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numerical technique is used to calculate the flow for each closure. However, after numerous 
attempts, convergent solutions for LT are not possible under the uniformly applied convergency 
criterion of normalized residual value less than 1CH for all calculations. On the other hand, if the 
convergency criterion is relaxed, convergent solutions become possible. Since all comparisons are 
made on the same basis, the results of LT are not listed in Table 4. 

Based on the near- wall analyses and the results listed in Table 3, it could be concluded that 
HL, PE and KLY could not perform well in their predictions of the limiting values at the wall 
because the 1 1, 22, 33 and 13 components of (D v + <J>P + O;: - £;;) are not zero at the wall unless 

ij U J J 

Ci is set equal to one or the random functions bj, and cj are required to satisfy certain 

constraints. Indeed, this is the case as indicated by the calculated limiting values shown in Table 4 

and the results are at variance with simulation data. First of all, the calculated normal stresses are 

fairly isotropic near the wall. This is evidence by the approximately equal values of a u , a v and a w 

calculated by HL, PE and KLY. The calculated a uv is several times larger than the data shown 

while the calculated a k is several times smaller. Even though the limiting values for (a 2 + a 2 + 

a w )/a k are approximately correct, the calculated values for k + /e + y+ 2 are not. This means that the 

closures are not internally consistent and asymptotically incorrect. On the other hand, S Y goes to 

zero at the wall for every component of (D^ + <D?. + d>:: - £::). However, this is achieved at the 

ij y J J 

expense of solving a modified form of Chien's £-equation 32 that does not satisfy the true boundary 
condition at the wall. Consequently, the limiting values calculated using this closure are all 
incorrect. This means that the £-equation of Chien 32 is not at all suitable for near-wall turbulent 
flow calculations. 

The other three closures, SH, LSH and LS, give very reasonable results compared to 
simulation data as well as measurements. This is a consequence of the fact that they are formulated 
to satisfy the asymptotic near-wall behavior of either the velocity-pressure-gradient correlation or 
the exact Reynolds-stress equations or both. The limiting values for (a 2 + a 2 + a^,)/a k are 
recovered exactly while the values for k + /£ + y +2 are correct to within 4%. There are two major 
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differences between these three closures though. One is in the predictions of a u and a k , while the 
other is in the calculations of a v and a uv . LSH gives a more correct predictions of a u and a k . 
However, its calculations of a v and a uv are greater than the data shown by approximately one order 
of magnitude. On the other hand, the reverse is true for SH and LS. There is one exception 
though, and that is the prediction of a v by SH. Its value is as large as that calculated by LSH. The 
reason could be due to the fact that LS tries to match the asymptotic behavior to higher orders of y 
for the 12, 23 and 22 components of (D^ + + Ojj - £jj). If the predictions of a u and a k are to be 

improved, one suggestion could be to try to match the 1 1 component of (D^ + OP + 0 V) - Ey) to y 2 
in the near-wall region. 

Another test of the internal consistency of the near-wall closures can be earned out by 
comparing the limiting values of the structure parameters, a^/a k , a^/a k , a^/a k and a uv /a u a v , with 
data. These values are also listed in Table 4 for comparison. Both the direct simulation results and 
measurements give an aj/a k value that is closed to zero at the wall. On the other hand, the value of 
aj/a k is about 3.5 to 6 times that of a> k and a uv /a u a v = 0.22 to 0.26, respectively. The measured 
ratio aj/a k is on the high side because of a probable incorrect measurement of a u . None of the 
closures, HL, PE, KLY and S Y, examined give limiting behavior of a^/a k , a v /a k and a w /a k that is 
even closed to the data, especially the calculated a^/a k , which is at least two orders of magnitude 
greater than the data shown. However, HL and PE give fairly accurate predictions of a uv /a u a v . 
The calculated ratio of aj/a k over aj/a k from SH is approximately equal to the measured result but 
its calculated a^/a k is two orders of magnitude larger than the simulated value and its prediction of 
a uv /a u a v is abouit half that of the measurement. As for LSH, its predictions of the ratio aj/a k over 
a^/a k and a uv /a u a v are twice as large as the data. On the other hand, LS's prediction of a v /a k is 
much closer to that of the data but its calculation of a uv /a u a v is about the same as SH. It seems that 
SH and LS are the only closures that could replicate the limiting behavior approximately. 
Therefore, it could be concluded that, overall, three closures are capable of reproducing the near- 
wall behavior more correctly than other closures examined. Furthermore, they are internally more 
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consistent in terms of the predicted limiting values of the structure parameters, the Reynolds 
stresses and the dissipation rate. 

The plots of -uv 7 , k + , and the three normal stresses u* , v* and w + are shown in 

rmi> tin s fills 

Figs. 1-3. In each plot, part (a) shows the near- wall distributions given in terms of the wall 
variable y + , while part (b) displays the overall profiles in terms of the outer variable (1 - y/H). 
Each panel of each plot shows the comparison of the modelled calculation with direct simulation 
result 2 ^. The calculations are represented by solid curves, while the simulation results are denoted 
by dashed curves. 

From these plots, it can be seen that LS gives the best overall prediction of ^uv+ in both the 
near-wall and outer region. The maximum -uv + is also calculated correctly (Fig. 1). The next best 
predictions are those given by LSH, KLY and HL. Here, the closures under-predict the maximum 
■uv + - The other closures, PE, SY and SH, on the other hand, give an incorrect prediction of the 
near-wall distribution of -uv + . In a fully-developed channel flow, the mean velocity only depends 
on the shear stress distribution across the channel. Therefore, if a closure under- or over-predicts 
the turbulent shear stress, its calculation of the mean velocity profile is more than likely to be 
incorrect also. 

As far as k + is concerned, HL, PE and SY under-predict its value all across the channel, in 
the near-wall as well as in the outer region (Fig. 2). HL and PE even gready under-predict the 
value of k + in the channel mid-plane. On the other hand, KLY over-predicts k + in the outer region 
but its prediction in the near-wall region is better than those of HL, PE and SY. The calculation of 
SH is, in general, slightly lowered than the simulation data. On the other hand, the near-wall 
prediction of k + by LSH agrees well with simulation data. However, it over-estimates the value of 
k + slightly in the outer region. Even though LS under-predicts the maximum k + just like LSH, it 
does give a fairly correct prediction of the k + profile across the channel. In general, the best 
overall predictions of k + are given by LSH and LS. 
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Finally, the same can also be said of the predictions of u^ s , vj^ and (Fig. 3). Only 
SH, LSH and LS could reproduce the anisotropic behavior of the normal stresses near the wall, the 
other closures invariably give a fairly isotropic prediction. Even though the normal stress 
predictions are fairly anisotropic for SH, its calculation of is substantially higher than data in 
the region y+ < 30. This over-prediction of v^, would lead to an incorrect calculation of the - uv 
budget and, thus, renders SH less attractive compared to LSH and LS. The major reason for the 
noted discrepancy lies in the modeling of the (d>y + Ojj) term. SH, LSH and LS attempt to 
partially model OP in the near-wall region. In the process, they also try to partially compensate for 
the incorrect near- wall behavior of Ojj. Consequently, all components of (D^j + OPj + Ojj - £jj) 
vanish at the wall independent of the behavior of the random functions b, and q.. This, 
however, is not true for HL, PE and KLY; therefore, it is not surprising that they fail to give the 
correct predictions of the normal stresses. Again, LSH slightly over-predicts the normal stresses 
in the outer region. 

The plots ofeH/uJ versus (1.0 - y/H) and eH/uJ versus y+ are given in Fig. 4. Here, the 
comparisons are carried out to assess the performance of the different closures relative to each 
other. With the exception of SY, the general trend of the predicted £ is about the same. The 
behavior can be described by a general increase of £ towards the wall. Just before the wall, a 
maximum is reached and this is followed by a rapid drop to a finite wall value. This trend, in 
particular, the near-wall behavior, is contrary to simulation data 26 , where £ increases to a 
maximum at the wall. Since all these closures satisfy the boundary condition that the wall £+ is 
equal to 2a k , the calculated behavior of £ in the near-wall region could not be due to an incorrect 
specification of the boundary condition. Rather, it is the consequence of inappropriate near- wall 
models for either the Reynolds-stress equations or the E-equation. However, judging from other 
predicted quantities, this incorrect behavior of £ in the near-wall region does not seem to have a 
significant effect on the calculations of the Reynolds stresses and structure parameters, especially 
in the case of LS, LSH and SH. 
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From the above comparisons, it is clear that HL, PE, KLY and SY perform poorly as near- 
wall closures. Therefore, their predictions of the Reynolds-stress and k budgets would also be 
incorrect. On the other hand, the budgets deduced from SH and LSH with the exception of the - 
uv budget are very similar while those of LS have been presented in Ref. 15. In view of this, 
only the budget plots of LSH and HL for u^ , k and - uv are shown for comparison with 
simulation data in Figs. 5-7. The sample plots presented here together with those given in Ref. 
15 serve to illustrate the discrepancies between closure calculations and data. It is hoped that, 
through this evaluation, an approach to improve the closures could be identified. 

All budget plots shown in Figs. 5-7 and those presented in Ref. 15 give an incorrect 
prediction of e^, and e, especially in the region, 0 < y + 20. According to simulation results, 
the absolute maximum of e j \ and e occurs at the wall. The model calculations, on the other hand, 
show that their absolute maximum appears away from the wall. There is, of course, a drop from 
these maxima to some finite wall values. Since viscous dissipation is balanced exactly by viscous 
diffusion at the wall, it follows that the calculated wall values of the viscous diffusion of u^ and k 
are also incorrect, in spite of the fact that their predicted trends are similar to those of simulation 
data (Figs. 5 and 6). In LSH's calculation of the budget, viscous diffusion is balanced by both 
dissipation and velocity-pressure-gradient correlation (Fig. 5). This is a consequence of the 
incorrect modeling of the velocity-pressure-gradient term. In general, all closures give a poor 
prediction of the behavior of the velocity-pressure-gradient correlation term in the near-wall region; 
in particular, the behavior in the 12 component (Fig. 7). This is especially true of HL where the 
calculated value is about four times larger than data. The reason is because v^ is responsible for 
production in the - uv equation and HL gives a substantially higher prediction of this nomal stress 
in the near-wall region (see Fig. 3a). Since production of - uv has to be balanced by the velocity- 
pressure-gradient correlation of - uv , the result is a totally incorrect prediction of this budget. 
Consequently, the - uv budget obtained from HL is shown in Fig. 7 with a scale different from 
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that of LSH and simulation data. This is another indication that HL is not an asymptotically correct 
near-wall closure. 

The wall values of e 12 , £22 and e 33 are zero and hence the viscous diffusion is also zero for 
these components. Even though these calculated trends are consistent with simulation data, the 
magnitudes of the calculated distributions of the viscous diffusion along the 12, 22 and 33 
components are several times larger than the simulated values. This evaluation, therefore, reveals 
that although SH, LSH and LS give good predictions of the turbulence statistics, their calculations 
of the Reynolds-stress and k budgets are still not quite correct. 

There are many reasons for this incorrect prediction. However, the two most important 
factors could be the e-equation and the modeling of the velocity-pressure-gradient term in the near- 
wall region. Lai and So^ did a preliminary investigation to show that the correct behavior of £j \ 
and e could be achieved by modifying the e-equation. Unfortunately, the modification also affects 
the calculation of other statistical properties. As a result, the overall predictions are not as 
satisfactory as before. In spite of this rather inconclusive result, their study did manage to show 
that, perhaps, the e-equation could be properly modified to remedy the incorrect prediction of £u 
and e. This approach has been adopted by So et al. 56 and they manage to successfully modify the 
e-equation to give a correct prediction of the turbulence statistics and the behavior of e in a flat plate 
boundary-layer flow. Further comments about the e-equation of So et al. 56 and its performance in 
terms of the prediction of the channel flow problem are given in the next section. 
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7 . Proposed Improvement 


The above analysis reveals that, among the eight closures examined, three manage to give 
the best overall predictions of the turbulence statistics in a fully-developed channel flow. 
However, not a single closure can reproduce the near- wall distributions of e^ and e correctly. 
Even SH, LSH and LS can only correctly replicate Ej j and e up to y + => 20. Thereafter, their 
predictions of and e are totally incorrect. Lai and So 15 argued that the inability of their closure 
to predict the near-wall distributions of e^ and e lies in the incorrect modeling of the near- wall 
production and dissipation terms proposed for the e-equation as well as in the near-wall modeling 
of the velocity-pressure-gradient correlation. Since the modeling of the e-equation is rather ad hoc, 
a first attempt to improve these predictions could be achieved by modifying the near-wall 
production and dissipation terms proposed for the e-equation. A preliminary effort has been made 
by So et al. 56 They used the e-equation of Ref. 15 as a base and suggested an alternative £, 
function based on the studies of Shima 12 and Mansour et al. 55 Hie proposed £ is given by 



In modifying the e-equation, So et al. 56 maintained the values of most model constants specified in 
LS but found that they have to change C E i to 1.50 and the definition of a to [1 - 1.5exp(-Re/lCH)]. 
Thus modified, the e-equation is incorporated into a k-e closure to calculate flat plate boundary- 
layer flows. Their results show that the calculated properties are in very good agreement with 
measurements and that the limiting values of the structure parameters and the near-wall e 
distribution behave similarly to the direct simulation result of a flat plate boundary layer 27 . 

In order to verify that the primary cause for the incorrect e^ and e behavior near a wall is 
the e-equation, it is suggested to improve one of the near-wall closures identified above by 
adopting the modified e-equation of So et al. 56 Such an improved near-wall closure could be 
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obtained from the Reynolds- stress models of Refs. 12, 14 or 15 and the e-equation of Ref. 56. 
Here, the example is carried out with the models of Ref. 15. The resultant closure, designated as 
LSZS, is used to calculate the fully-developed channel flow considered above and the calculated 
limiting values of the Reynolds stresses and structure parameters are listed in Table 4 for 
comparison. Furthermore, the plots of the Reynolds stresses, k and £ are shown in Figs. 1 4. 

Near-wall budgets of -w and k are plotted in Figs. 5 - 7. It can be seen that LSZS gives the 
best overall prediction of the limiting values of the Reynolds stresses and structure parameters (see 
Table 4). Most important of all, the calculated a^ agrees to within 4% of the simulated data. This 
means that, for the first time, the wall e + is predicted correctly. In addition, the calculated a u and 
a w are in very good agreement with data compared to other closures, including LSH. As a result, 
the calculated structure parameters are also in better agreement with data. 

A comparison of the distributions of the Reynolds stresses, k and e reveals that LSZS gives 
better predictions of the overall profiles, including the near-wall behavior of e (Figs. 1-4). In 
spite of the good agreement shown, two slight discrepancies should be noted. One is in the near- 
wall prediction of the distribution of -"uv, while another is in the calculation of k near the channel 
center. LSZS over-predicts -uv in the near-wall region and is a conseqence of the incorrect 
prediction of a uv , which is probably due to an incomplete modeling of the pressure diffusion term 
near the wall. It also under-predicts k near the channel center and is primarily caused by an under- 
estimation of v 7 in the outer region of the channel. This under-estimation, however, is most likely 
tied to the model for pressure redistribution. Therefore, the calculated results could be further 
improved by modifying these two models. 

For the first time, the near- wall behavior of £n and e is predicted correctly in a qualitative 
sense (Figs. 5 and 6). The calculated distributions do not show a distinct plateau in the region, 8 < 
y+ < 15 , compared to those given by direct simulation. However, they do increase to a maximum 
at the wall. The calculated wall values of -e n and -e are -0.28 and -0.172, respectively, and they 
compare favorably with the simulated values of approximately -0.26 and -0.166. Consequently, 
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the wall values of viscous diffusion for the u 2 and k components are also predicted correctly, and 
so is the steep rise of the viscous diffusion term in the near-wall region. Even though the 
modification improves the prediction of the budgets of u 2 " and k, it has very little effect on the -uv 
budget (Fig. 7 ). Again, the calculated -E^ is 2 to 3 times larger than the simulated value and is in 
close agreement with that predicted by LS. As for the budgets of v^ and w 2 , they are essentially 
identical to those of LS and, therefore, the predictions of £22 and £33 are again several times larger 
than the simulated values. In view of this, the improvement fails to completely remedy all the 
discrepancies noted in Section VI. 

From the above comparisons, it can be seen that the improvement made to the £-equation 
does not significantly alter the calculated turbulence statistics (Figs. 1 - 4 ) and the various terms in 
the Reynolds-stress and k budgets (Figs. 5 - 7 ). It only changes the distributions of e n and £ near 
the wall (Figs. 5 and 6). This means that the behavior of and e is very much influenced by the 
modelled £-equation and substantiates the approach taken by So et al. 56 to improve near-wall 
closures. Therefore, this analysis suggests that further improvement to second-order near-wall 
turbulence closures is still needed and this could be achieved by modifying the models proposed 
for the velocity-pressure-gradient correlation term. It is hoped that, through this further 
improvement, the behavior of £j2, £22 and £33, the near-wall distribution of -"uv and the centerline 
value of k could be predicted satisfactorily. 

Based on the above improvement, it is now possible to calculate the turbulence statistics of 
a fully-developed channel flow fairly correctly using a second-order near-wall turbulence closure. 
Also, for the first time, the behavior of Eu and £ in the near-wall region has been calculated 
correctly. However, improvements are still needed to remedy the discrepancies noted in the 
predictions of E^, £22. £33 and -Irv in the near-wall region and the calculation of k near the 
centerline. 
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8 . Concluding Remarks 


A total of eight near-wall Reynolds- stress closures has been reviewed. All closures 
examined essentially adopt the same high-Reynolds-number models for turbulent diffusion, 
dissipation and pressure redistribution. Invariably the pressure diffusion part of the velocity- 
pressure- gradient correlation is neglected in all high-Reynolds-number closures. Different 
approaches have been proposed to modify these closures for near-wall flows. Some adopt the 
philosophy that only the model for dissipation needs modification, while others argue for 
modifications of the models for both dissipation and pressure redistribution. Some others also 
argue for the inclusion of pressure diffusion effects in the near-wall region; therefore, suggesting 
that the velocity-pressure-gradient correlation has to be modelled in the near-wall region rather than 
pressure redistribution alone. All modifications are proposed so that the resultant closure satisfied 
the constraints imposed by the near-wall asymptotic behavior of the Reynolds stresses. With the 
exception of SY, all closures examined solved a transport equation for the true dissipation rate of 
the turbulent kinetic energy. On the other hand, S Y solves a pseudo dissipation-rate equation that 
satisfies a homogeneous wall boundary condition. These dissipation-rate equations are further 
modified to mimick the exact near-wall behavior. As a result, a uniform claim of all these closures 
is that they are asymptotically correct as a wall is approached. 

The modelled equations are examined for their asymptotic behavior near a wall and the 
results are compared with the behavior deduced from the exact Reynolds-stress equations. Among 
the closures examined, it is found that four of them (HL, PE, KLY and LT) fail to go to zero at the 
wall for the 11, 22, 33 and 13 components of the equations without imposing some kind of 
constraints on Ci or on the behavior of the random functions aj, bj and q. The others approach 
zero at the wall independent of the behavior of the random functions; however, their asymptotic 
behavior is not consistent with that required by the exact equations. LS matches the behavior to a 
higher order of y than SY, SH and LSH for the 12, 23 and 22 components (see Table 3), while the 
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behavior of the 1 1, 33 and 13 components is similar to that of SY, SH and LSH. A further check 
on the internal consistency of the closures is carried out by using them to calculate a fully- 
developed flow in a two-dimensional channel at a Re = 6500. These results are then compared 
with direct simulation and experimental data 25 ’ 26 * 45 . Specifically, the limiting values of the 
Reynolds stresses, the structure parameters and the dissipation rate are assessed. A common 
numerical technique is used to perform the calculation for each closure and a uniform convergency 
criterion is imposed on all calculations. Under these conditions, LT fails to give convergent 
solutions to the governing equations. Of the remaining seven calculations, it is found that four 
closures (HL, PE, KLY and SY) fail to give limiting values that approximate the simulation and 
experimental data examined. Their failure is most pronounced in the prediction of the turbulent 
stress normal to the wall and the shear stress. For example, the data examined shows that a^/a k 
goes to zero at the wall, while the closure calculations of HL, PE, KLY and S Y give finite values 
that are comparable to the other structure parameters at the wall. On the other hand, SH, LSH and 
LS give internally consistent results that are asymptotically more correct compared to other closure 
calculations. However, one drawback of LSH is its calculation of a^/a k , which is still one order of 
magnitude larger than that obtained from LS, while the values of a u and a k deduced from SH and 
LS are quite a bit less than those of the simulated values. Therefore, SH, LSH and LS have their 
defects and should be further improved. 

A comparison of the closure calculations with data further reveal that, with the exception of 
LSH and LS, all closures examined either over- or under-predict the turbulent shear stress in the 
near-wall region or under-predict the maximum shear. This means that they are more than likely to 
give an incorrect calculation of the mean velocity for the channel flow considered. Even though LS 
under-predicts the maximum turbulent kinetic energy, its overall prediction of the k profile is in 
good agreement with data. On the other hand, SH under-predicts k across the whole channel while 
LSH over-predicts k in the outer region. The other closures essentially give an incorrect prediction 
of k across the channel. Finally, only SH, LSH and LS are capable of reproducing the anisotropic 
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behavior of the normal stresses near the wall, while the other closures give a fairly isotropic near- 
wall behavior. Even then, SH's prediction of v 2 in the near-wall region is about twice as large as 
the simulation data. Consequently, its calculation of the -uv budget is greatly in error. The 
reason is traced to the imbalance of the modelled equations in the near-wall region. In spite of 
these good comparisons by SH, LSH and LS, their calculations of e n and e in the near-wall 
region are totally incorrect. The three closures give a maximum e n and e away from the wall and a 
rapid drop to a finite value at the wall. On the other hand, simulation data shows that e n and e 
reach their maximum value at the wall. As a result, the wall values of En, e and the viscous 
diffusion of u 2 and k are greatly under-predicted. 

A conjecture is then made that the incorrect prediction of e n and e could be remedied by 
improving the near-wall behavior of the e-equation. This conjecture is essentially verified by the 
calculations of LSZS . The LSZS closure is made up of the Reynolds-stress models of Ref. 15 
and the e-equation of Ref. 56. With this improvement, the calculated results are in much better 
agreement with data, including the near- wall distributions of e t j and e. The calculated limiting 
values of the Reynolds stresses give the best overall agreement with data, and so are the structure 
parameters. Also, for the first time, the wall values of e^ and e are replicated by LSZS. This 
means that the wall values of the viscous diffusion of u 2 and k are also calculated correctly. 
However, some slight discrepancies exist in the predictions of e^, £22’ e 33 an< * " uv * n the near ' 
wall region and k near the channel center. These discrepancies arc relatively less important in the 
overall comparison. They could be attributed to an incorrect modeling of the velocity-pressure- 
gradient correlation near the wall and pressure redistribution near the channel center. Therefore, 
further improvements are required if a complete replication of the simulation data by a statistical 
model of turbulence is desirable. 
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Table 2. Summary of near-wall turbulence properties. 
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Based on boundary layer momentum thickness 
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Figure la. Distributions of -uv + in a fully -developed channel flow: model calculation 

direct simulation data 25 ; near-wall behavior. 
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Figure 2a. Distributions of k + in a fully-developed channel flow: model calculation 

direct simulation data 25 ; near-wall behavior. 


48 


o 

in 


o 

LO 


0.0 0.2 0.4 0.6 0.8 1.0 

1-y/H 


«r LSH 


o 

+ « 

o 

c\i 



0 0 0.2 0.4 0.6 0.8 1.0 

1-y/H 

■ i KLY 



0 0 0.2 0.4 0.6 0.8 1.0 

1-y/H 

5 i HL 



0 0 0.2 0.4 0.6 0.8 1.0 

1-y/H 


° 0.0 0.2 0.4 0.6 0.8 1.0 

1-y/H 

S 1 SH 


0 0 0.2 0.4 0.6 0.8 1.0 

1-y/H 

h SY 


0 0 0.2 0.4 0.6 0.8 1.0 

1-y/H 

S i PE 


0 0 0.2 0.4 0.6 0.8 1.0 

1-y/H 


Figure 2b. Distributions of k + in a fully-developed channel flow: model calculation 

direct simulation data 25 ; outer flow behavior. 
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Figure 3a. Distributions of uj t ’ ms , v^ ns and w^ s in a fully-developed channel flow: model 

calculation; — , direct simulation data 25 ; near-wall behavior. 
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